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Gauge invariance of scalar perturbations is studied together with the associated equations of 
motion. Extending methods developed in the framework of hamiltonian general relativity, the 
Hamilton-Jacobi equations are investigated into the details in Loop Quantum Cosmology. The 
gauge-invariant observables are built and their equations of motions are reviewed both in Hamil- 
tonian and Lagrangian approaches. This method is applied to scalar perturbations with either 
holonomy or inverse- volume corrections. 
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I. INTRODUCTION 

Loop Quantum Gravity (LQG) is a non-pcrturbative 
and background-independent quantization of General 
Relativity (GR) (see [l[ for reviews). Recently, it has 
been reahzed that different vievifs, based on canonical 
quantization of GR, on covariant quantization of GR and 
on formal quantization of geometry lead to the very same 
LQG framework. Although other approaches are still 
much debated, this makes LQG a very promising model 
to address the outstanding question of quantum gravity. 

The application of LQG ideas to the universe as a 
whole is called Loop Quantum Cosmology (LQG) (see 
for reviews). This is basically the symmetry- reduced 
version of the theory. So far, LQG proved to be interest- 
ing both as a model of the early universe, solving the Big 
Bang singularity, and as a way of possibly testing LQG 
ideas. At the effective level, LQG modifies the usual 
paradigm by two main corrections: the inverse-volume 
terms, basically arising for inverse powers of the densi- 
tized triad, which when quantized become an operator 
with zero in its discrete spectrum thus lacking a direct 
inverse, and holonomy corrections coming from the fact 
that loop quantization is based on holonomics, rather 
than direct connection components. 

To investigate the observational consequences of 
those LQC-induced modifications, it is most useful to 
construct rigorously gauge-invariant variables. It is well 
known, even in standard GR, that among the solutions 
of field equations for perturbed variables, some are 
unphysical modes corresponding to a mere coordinate 
transformations. 

In this article, we basically extend the method intro- 
duced in 0]. We start with first order constraints in 
the FLRW metric. Then, using a generating function, 
the variables (perturbations of the densitized triads 
and their conjugate momenta) are changed according 
to {SK,5E) — > {jrm'^m)- The first order constraints 
are re-expressed with respect to (7m, tt^). The gauge- 
invariant variables {Q,P) are obtained thanks to a 
natural generating function S and the dynamics is 
derived through anomaly-free second order constraints 



in terms of {Q,P). Then, the Mukhanov variables v, R 
and z are given. Finally, the method is applied to the 
case of LQG with both holonomy and inverse-volume 
corrections. 

This approach exhibits several advantages: 

• the treatment is purely Hamiltonian with easy com- 
putations, 

• the Mukhanov variables v and R arc obtained di- 
rectly and the equation of motion is easily found 
without using Bardeen Potentials, 

• it helps to construct an anomaly-free algebra by 
imposing relations on the Poisson brackets, 

• the z variable can be found without ambiguity and 
in a quite simple way, 

• the generating functions are clearly defined, easy 
to handle and allow one to trace back deeply the 
origin of gauge invariance, 

• it works for any kind of constraint theory. 

The paper is organized as follows. In the two first 
sections, we introduce the framework of Loop Quantum 
Gosmology and some elements of analytical mechanics 
useful to implement the Hamilton-Jacobi method. Then, 
we show the main steps of the proposed procedure and its 
application to the cases of holonomy and inverse- volume 
corrections. 



II. LOOP QUATUM COSMOLOGY 
FRAMEWORK 

In General Relativity, when the ADM formalism is cho- 
sen, space-time is foliated into a family of spacelike 3- 
surfaces and the dynamics is given by constraints. The 
fundamental variables are the space metric qab, together 
with N, the lapse function, and N"", the shift vector, 
which describe how each of the "leaves" of the foliation 
are welded together. The metric is written as: 



-N^dt^ + Qabidx" + N''dt){dx'' + N^dt). (1) 
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In the LQC formahsm, the spatial metric is expressed in 
terms of triads that arc related to the spatial metric 
by: 



Qab — ^a^b- 



(2) 



The first basic variable (for a detailed introduction, see 
[l[) is the Ashtckar connection: 



(3) 



where 7 is the Barbero-Immirzi parameter, T\ is the spin 
connection and /C* is the extrinsic curvature. The second 
one is the densitized triad: 



Et = [det el) e^ 



(4) 



The conjugate variables follow the symplcctic structure 

{Al{x),E'^{y)}^K^S]5y{x-y), (5) 

where n = 8t:Q. The canonical Einstein- Hilbert action 
in this formalism reads as 



S 



EH 



dt 



(fa 



A^E^ - G[A''] - D[N°] - H[N] 



(6) 

where G'[A*] is the Gauss constraint, D[N"'] is the 
diffcomorphism constraint, and H[N] is the Hamiltonian 
constraint. The diffcomorphism constraint generates 
deformations of a spatial slice so that, when it is satis- 
fied, spatial geometry does not depend on the choice of 
space coordinates. General covariance of the spacetime 
geometry (including the time coordinate) is ensured by 
the Hamiltonian constraint. Finally, as a set of triad 
vectors can be rotated without changing the metric, 
there is an additional SO (3) gauge freedom. Invariance 
of the theory under those rotations is guaranteed by the 
Gauss constraint. This latter constraint will be solved 
explicitly by the parametrization we use for the variables. 

Taking into account perturbations in a FLRW uni- 
verse, one has to deal with the perturbed spatial metric 
S"tab such that: 



qab = a'^iSab + hab), 



(7) 



where a{t) is the scale factor. The perturbed Ashtekar 
variables will then be related to the perturbed metric 
and it is straightforward to see that the background and 
perturbed densitized triad obeying Ef ~ Ef + SEf are: 



and 



SE:^-p{-5jt + dj^-S1) 



(8) 



(9) 



On the other hand, the extrinsic curvature is given 

by 



(10) 



The homogeneous and isotropic background (p, k) satis- 
fies 

K 



3' 



(11) 



and the perturbed part {5Ef,SKj) fulfills 

{6Kl{x), SE'^iy)} = KS]S'j'ix - y). (12) 

In 0], the original variables are the spatial perturbed 
metric S'^ab and its conjugate momentum Stt"-^. As SEf 
is linear in 5^, as it can be seen in Eq. ([S]), it is possible 
to follow the same procedure, with only minor modifica- 
tions due to the fact that now the fondamental variables 
arc SE and SK. 

The study of the homogeneous and isotropic universe is 
an important first step for any tentative theory of quan- 
tum cosmology. In the framework of LQC, this led to the 
famous replacement of the Big Bang by a Big Bounce. 
Investigating perturbations is the next logical step to 
probe possible deviations from the standard model. This 
has already been studied in many articles (see, e.g., Q), 
especially for gravitational waves and subsequent con- 
sequences of the B-mode spectrum of the Cosmological 
Microwave Background (CMB). 

We now turn to the study of scalar perturbations of the 
metric ((/>, ip, B and E) that arc observationally relevant 
as they can be used to compute the well measured tem- 
perature CMB spectrum. The perturbed FLRW metric 
in conformal time can be written as 

d.s^ = 0^(77) [-(1 + 24>)drf + 2daB ■ dx^drj 

+ ((1 - 2iP)6ab + 2dadbE)dx''dx''] . (13) 

Comparing this expression with Eq. ([T]), one obtains the 
perturbed lapse function and perturbed shift vector as: 

SN = N(j>; SN" = d^B. (14) 

Using the definition of the densitized triad, one can also 
see that 



SEf ^ p i^2^Sf + {SI A - d,d'')E) . 



(15) 



Starting from Eq. (|15p . all the useful equations will be 
derived using a clear algorithm. 



III. HAMILTON-JACOBI EQUATION 

This section is heavily based on [5|. 

A. Canonical transformations 

When dealing with general transformations of coordi- 
nates, one has to consider the simultaneous transforma- 
tions of independent coordinates and momenta, qi and pi 
to a new set Qi and Pi, through (invertible) equations : 



Q^ 
P, 



Qiiq,p,t), 

Pi{q,P,t). 



(16) 
(17) 
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Theses equations define a transformation from a point in 
phase space to another one. 

In the Hamiltonian mechanics framework, only those 
transformations for which the new Q,P are canonical co- 
ordinates are interesting. This means that there exists a 
fimction K{Q,P,t) such that the equations of motion in 
the new set are in the Hamiltonian form: 



dK 



p, 



dK 



(18) 



Transformations for which Eqs. ([T5|) are fulfilled are said 
to be canonical. 

The function K plays the role of the Hamiltonian in the 
new coordinate set. For the treatment to be fully generic, 
for all systems with the same number of degrees of free- 
dom, Eqs. (HHI) must be the equations of motion in the 
new coordinates and momenta whatever the initial form 
of H . No matter wether one deals with an harmonic os- 
cillator or with a two-dimensional Keplerian problem. 
If Qi and Pi are to be canonical coordinates, they must 
satisfy the "modified" Hamilton principle 



(19) 



5 j (P^Q^- KiQ,P,t))dt = 0, 
whereas, as usual, 

6 [ {p,q^- H{q,p,t))dt = 0. 



(20) 



Both equations will be satisfied if the integrands are con- 
nected by the relation: 

dF 

X{M,^H{q,p,t))^P,Q,-K{Q,P,t) + —, (21) 



and both sets verify the Poisson Bracket: 
{q,p} = {Q,P}^l. 



(22) 



F is useful when mixing half of the old variables with 
the new variables and will then be considered as a bridge 
between the two sets of canonical variables. It is called 
the generating function of the transformation. On can 
define 4 such generating functions: 



Fi{q,Q,t), 
F3{p,Q,t), 



F2{q,P,t), 
Fi{p,P,t), 



with the following properties: 



dFi 
dqi 



Pr = - 



dF\ 
dQi 



dFi 
dt - 



K = H 



m 

dqi ' 



dP,: 



dt - 



(23) 

(24) 
(25) 

(26) 
(27) 



It 
K 



dpi ' 



= H 



dt ' 



dF4_ 
dpi ' 



dF4_ 
dP' 



K ^ H 



dFi 
dt 



(28) 
(29) 

(30) 
(31) 

(32) 
(33) 



In the particular case where, for instance, 
F3{p,Q,t) = -p-Q, 

Eq. (UHl) gives 

g = 0, P = p, 

which corresponds to the identity transformation. 

B. Hamilton- Jacobi equation 



In order to solve a problem of mechanics, it is useful 
to formulate it with the best suited variables, for ex- 
ample the angle-action variables. Then, one solves the 
Hamilton- Jacobi equation written thanks to a generat- 
ing function, S, which changes initial variables to new 
appropriate coordinates. The equation is basically given 

by 



H qt,Pi 



dS 



Ui 



(34) 



For gravity, one has to deal with constraints, as intro- 
duced in the previous sections. General methods to solve 
the Hamilton- Jacobi equation in this case, with = 0, 
arc given in Q. 



IV. FIRST CHANGE OF VARIABLES (7„,7r„) 

In this section, we show in some details the way to 
proceed in order to find easily the gauge-invariant quan- 
tities. This "algorithm" of resolution, originally used in 
Q but not fully detailed, can be applied in many situa- 
tions where perturbations are considered (see for instance 
3). The Hamilton- Jacobi framework has already been 
extensively studied and used in General Relativity, as, 
e.g.. in but the method presented here focuses on the 
goal of directly deriving some gauge-invariant variables 
useful for observations. 

It will here be explains for the case of General Relativ- 
ity but expressed with variables that can be further used 
in the framework of Loop Quantum Gravity, as investi- 
gated in the last section of this article. 
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A. new variables 

Following 0, we will define, from an appropriate gen- 
erating function, two "new" variables 7 and their con- 
jugate momenta tt, related to SK and SE, so that the 
equations get simplified. This is nothing else than re- 
expressing the perturbations E and ip. Fundamentally, 
this does not bring any new physical information as it is 
possible to obtain the very same results starting directly 
from -0 and E . 

We will use the Fourier transformed variables such that, 
for instance. 



SE^{k,t) ^ / d^x e"^ '' 5E^{x,t), 



(35) 



leading to: 



5Et{k,t) = j d^xe-''''-'p (^-2^(6^) + {S^ - ^)k^E^ . 

(36) 

Working in the Fourier space greatly simplifies most 
equations and can add some freedom. In our case, we 
define two vectors A"* (m=l,2) in the Fourier space such 
that: 



A„,- 



so, 

ad, , 



k k'^ 



(37) 
(38) 



The variables a and b depend on the choice of the basis 
but, as we will show later, the final results do not depend 
on them. The scalar product of these vectors is propor- 
tional to 2fc2. This is in sharp contrast with the situation 
studied in Q where one had: 



Ai 



(39) 
(40) 



The difference is due to the choice of the perturbation: 
in our case we use E, whereas the "standard" /j, was used 
in leading to: 

6Enk,t) - J d^xe-^'-^p (^-2V^(<5n + (^ - ^'Jf)^^ 

It is easy to see that A\ ■ Aj^ = 0. However, this is not 
in principle necessary and both approaches arc strongly 
related and lead to the same results. Having defined 
theses vectors, instead of working with ip and E, we will 
use two other variables 7,„ (m = 1,2) such that: 



SEfik,t)=j,Al,+j2Al, 



(41) 



As suggested before, these new variables are just related 
to the perturbations through 



071 
&72 



2pV, 
pk^E. 



(42) 
(43) 



Using Eq. (HI]), one can express them in terms of 6E 
such that: 



«7i = '-^SEf, 
b^.^-lU'-fsEt-SE:^ 



(44) 
(45) 



Taking the trace of (|4T|) indeed leads to: 

SEj = 3a7i + 2672, (46) 

and expressing 71 as a function of 72 in Eq. (j4ip . one 
obtains: 



SEt = lsE',6t + 572 [Ist - ^ 



This can be expressed as: 



1 



SEl' = -SEjAi + bj2Al, 



(47) 



(48) 



by replacing S^ij (as used in 0]) by SEf. This is the first 
bridge between the two approaches. When solving this 
equation by multiplying by A\^^ , one naturally obtains 
Eqs. (gH) and 

Furthermore, wc can show that both approaches are in 
fact equivalent. Starting from one, for example using 
/i. wc can derive the expression of when E is used 
in terms of 7^. This is simply performed by using the 
relation Eq^^® which relates 5"fij to 5Ef. To be consis- 
tent with [3|, we re-define (due to our conventions) the 
variables as: 

p5^t =p^^Al+p^^Al. (49) 
Noticing that, from Eq. ([5]), 

pS^t = 5Ei5<} - 25Et, (50) 
we can re-express the approach of Q as: 



7i 



(51) 



' li^^^^"-^'^^^^)- (52) 



From Eqs. ©, (gT]), (gl]) and it follows that 



071 5f + 2by, 



Taking the trace gives: 



k^kfi 



26 



Pli = 0.11 + y72, 



and Eq. ([55)) becomes 



PI2 



k 3 ^ 



(53) 
(54) 

(55) 
(56) 
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This leads to the expected equations: 
2&72 



1 



071 = P{li - -^12 



(57) 
(58) 



Re-expressing them as functions of SE, with Eq. (lO; 
([5T|) and ^ leads to Eqs. (04]) and (|45l) . 



B. generating function S^, 

In the framework of perturbed LQC, the canonical 
variables are 5K and their conjugates 6E. We have seen 
that it is possible to make a transformation from SE to 
7m. To this aim, we have to define the corresponding 
conjugate variables tt^ of 7^, which will depend on 6K. 
As reminded above, there exists 4 generating functions 
allowing one to define new sets of variables. In our case, 
we define the momentum with a generating function S-y 
such that 



97m' 



(59) 



where S-y is then a function of As 7 is a function 
of SE, the momenta do not depend on SE and we can 
therefore set 



s-y = cskia:. 



\lr. 



(60) 



where c is a constant. At this stage, we might consider 
two cases. First, one may chose to have 7™ as canonical 
coordinates and iTm as their conjugate momenta: 



SK = q, SE ~ p, 



(61) 



In this case, will be similar to a /i((?, Q) function and 
the conjugate momenta are: 



dim 



-cskia: 



(62) 



In the second case, as 7™ are related to SE = p, one might 
want to have now tt^ as canonical coordinates and 7^ as 
their conjugate momenta: 



SK = q, 



SE^p, 



In 



(63) 



In this case, S-y will be similar to a /2(<7, P) function and 
the conjugate momenta are: 



dim 



cSKlA':. 



(64) 



Comparing both cases, one can see that changing c — > — c 
exchanges one case for the other one. From now on, we 
focus on the first case: the 7^ will be the canonical coor- 
dinates, and -Era their conjugate momenta. Nevertheless, 



considering either 7 = Q or 7 = P , the algorithm natu- 
rally leads to the same correct gauge-invariant variables. 
In our choice, those variable are precisely the Mukhanov 
variables. 

Remaining as general as possible, one can finally write: 



TTi = -acSKj 

71-2 = 



cb[SKi-^SKl 



(65) 
(66) 



It is useful for the following computations to re-express 
SK as a function of 7r„, such that 



where 



K = b,si + b2 



k^ka 
k^ka 



(67) 

(68) 
(69) 



Multiplying Eq. (|57)) by — cA™ leads to conditions on ai, 
&i, 02 and 1)2 through 



■Km = -CAZ (a>i + fo>2) , 

and, consequently. 



(70) 



It is useful to multiply the previous equation by SEf so 
as to obtain: 

SKlSEf ^ -^S^.i + ^fs^-^^V. 



2bc 



fc2 



1 



(7i7ri -I- 727r2) 



(72) 



As explained in the next section, this might be inter- 
preted as a conservation equation. 



C. Poisson brackets 

We now have defined the conjugate momentum tt™ of 
7m- As for the original Ashtekar variables, these new ones 
will obey some Poisson bracket relations. Going through 
the computation leads to: 



{7m7 7r„} = K c S„ 



(73) 



The transformation can be said to be canonical as the 
variables have a symplectic structure such that new Pois- 
son bracket is related to the old one through 



{l,n,TTn} = {SKi^{x),SEf{x)} C S„ 



(74) 



In the next sections, we will consider a universe filled with 
matter, and, in particular, with a massive scalar field ip, 
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with its conjugate momentum tt, so that {<^,7f} = 1. For 
simphcity and without any lack of generahty, we therefore 
set KC = 1. So, 



{jrn, 7r„} = 6r, 

Using Eq. ((72)) . we see that 



SKlSEf 



7m TT,, 



{7m, 7r,ri} 



(75) 



(76) 



As we have chosen the 7m, related to 5E, as the canoni- 
cal coordinates (instead of 6K in the usual theory), this 
leads to the appearance of a minus sign in the previous 
equation. 

The new set of symplectic variables, 7„i and tt„i, is now 
well defined as: 



bj2 = -TT 3 — 



SEf - SEfSl 



-acSKi 



TT2 = -bclSK^- 



SIC 



Their Poisson brackets are given by: 

{71,7^1} = {72,71-2} = 1, 
{71,7^2} = {72,71-1} = 0. 



(77) 
(78) 

(79) 
(80) 



(81) 
(82) 



As explained before, at this stage, nothing new emerges. 
This transformation is just useful to obtain simpler equa- 
tions. 



D. First order constraints in term of (7a,7ra) 

Gauge-invariant variables are derived from first order 
constraints. To use our new set of variables, we now 
have to re-express the constraints in terms of 7™ and tt„i. 
With the Ashtekar variables, for a universe filled with a 
massive scalar field (p (with conjugated momentum tt), 
the diffeomorphism constraints in the ADM formalism 
reads as 



DIN"] 



(fa 



(83) 



at first order in constraint densities. In fact, as in this 
case N"" = 0, only the 2?'^^ term remains. Its gravita- 
tional and matter components are: 



(1) 



1 



- {~kS%iSEt)+pdSKi)-pddiSK^)) 



2?« = nidjip). 



(84) 
(85) 



As far as the Hamiltonian constraints are concerned, 
one has: 



H[N] 



(86) 



where the first order constraint densities are: 



-1/(1) _ 

tin — 



2k 



VP 



-j=dcd'SE'^ 
VP 



515 E-T: 



ttSt: 



53/2 2p3/2 2p 



-3/2 



5i5E'i 



2p 



(87) 
(88) 
(89) 



As in g, we define 70 = Sip, and ttq ~ Stt, such that 7^ 
(a = 0, 1, 2) correspond to the old canonical coordinates 
"q" and {7a,7ra} = 1. The expressions of 7a represent 
the " maximal" and " fundamental" decomposition of the 
perturbations. What was done so far is nothing else than 
a decomposition of the theory in terms of those pertur- 
bations. 

After a Fourier transformation, and using Eqs. (j4ip and 
(j7ip . both first order constraints can now be expressed 
as functions of (7a,7ra) such that: 



H^^'>[5N]^5N (n'^a'' +n 



_ 5N 

~Tv 

+pV'7o 



2pk TTl ^ 

h Lp-KQ 

CK a 



ail 

K 



/(I) 
M 

K 



3P 



[2P 
3 



-Kip 



(90) 



and 



fc OTi ^ 

p K 



iPikaSN") 



^IPJQ 



1 TT2 
CK b 



(91) 



where V' and V" refer respectively to the first and second 
derivative with respect to the scalar field ip. The nota- 
tion "H'^' [5N] and 2?'^^ [5N'^] are used, in agreement with 
most papers, to emphasize that those expressions are in 
fact second order ones due to 5N and 5N°' factors. We 
have also simplified the results by using the Friedmann 
equation (calculations are derived in the Appendix): 



k' = 



pV 



(92) 



and the fact that the equation of motion for the back- 
ground variables reads as: 



ip : 



TT 
P' 



(93) 
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Wc have thus expressed the first order constraint densi- 
ties (PO)) and (|nT|) , as functions of the new set of symplec- 
tic coordinates (7a, tTq). In the next section, we will show 
that it is possible to make a last transformation toward a 
new set of coordinates {Q,P), meaningful for cosmology 
and gauge-invariant, using a generating function S and 
the Hamilton- Jacobi equations. 



V. GAUGE TRANSFORMATION AND THE 
MUKHANOV-SASAKI EQUATION 

To describe physical effects, one has to deal with gauge- 
invariant quantities. The goal of this section is to address 
the specific issue of gauge invariance within the canonical 
formalism. 



A. gauge-invariance in the canonical formalism 

In a canonical formulation, gauge transformations are 
directly generated by Poisson brackets of the fields with 
the constraints. In the covariant language, gauge trans- 
formations are studied as perturbation transformations 
under the coordinate change 

^x"' ^x" (94) 

generated by vector fields Under this coordinate 

transformation, any tensor field receives a correction 
given by its Lie derivative along As defined in Q, the 
part of the transformation relevant for the scalar modes 
can be parametrized by two scalar functions and ^ 
such that 



formations of basics variables as: 

Si^o^^^SEt = 2pH„St - PiS^k^ - kak')i, (100) 

+ f (-f +P^)?o<5:, (101) 
5i^o^^]6lp = if£,o, (102) 
%o^^]j7r = -p0k\ - p^V'^a- (103) 

With these expressions and the definition of 7,„ and 
TTm, it is easy to see that: 

- 0, 
= 0, 

which means that the first order constraints (|90p and 
([5T|) are gauge-invariant. Another way to see this is to 
compute directly the Poisson brackets: 

{if(2)[j7Vi],7?(2)[^iV2]} = 0, (104) 
{F(2)[5Ar],z?(2)[^iV^} = 0, (105) 
{d'-^^[SN^],D^'^'^[SN^]} = 0, (106) 

and replace the SN and JiV" by their equivalents. 
What is shown here has been noticed in Q. As it will 
be emphasized in the next section, this means that to 
obtain gauge- invariant quantities, the algebra should not 
only be anomaly-free, but should also have a null first 
order Poisson bracket. 



e' = (e°,a^O- 



(95) 



Along this vector, a variable X will undergo a transfor- 
mation given in the canonical formalism by: 



where 



(96) 

(97) 
(98) 



It is easy to relate the canonical approach to the Lie 
derivative by noticing that: 



{X + SX, D[C]} = Cf {X + SX). 



(99) 



In the framework of LQC, using the densitized con- 
straints (|84j|89p in Eq. one expresses the trans- 



B. gauge-invariance with the Hamilton- Jacobi 
equation 

In the Hamilton- Jacobi equation, the momentum is 
expressed in terms of a generating function S and a 
new transformation is performed. As stated in 0], there 
arc differences between the classical case where standard 
Hamiltonians are used, and the case studied here where 
we rely on constraints and reduce the phase space. In the 
latter case, the Hamilton- Jacobi-like equation has to be 
directly solved. As the total first order constraint (den- 
sity) has to be null for all SN and SN°-, which play similar 
(although slightly different) roles, one can separate the 
equations and solve the two Hamilton- Jacobi-like cxprcs- 



n'-'HSN] (^7„,7r„ = |^^ = 0, (107) 
2?(''m(^7a,T„ = |^^ = 0. (108) 
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Because 'H^'^\5N] and V'^'^'i[5N''] are linear in (7a,7ra) in 
((50)) and (pij) . the more "natural" generating function to 
consider is a quadratic function S ~ f2{la = <1-,P) such 
that: 

S = ^A„/37a70 +Sa7a-P2, (109) 

where A^p is a 3 x 3 matrix. Taking into account the 
properties of the generating function, the conjugate vari- 
able of P2 is given by: 

Q2=gJS-= Bala- (110) 

To show where the gauge-invariance of the new variables 
Q2 and P2 comes from, one can synthetically write that 

= {5N) {EaTTa + ^blb) 

= {5N){{EaAab + ^b)lb+EaBaP), (111) 

and 

= ^ • {k^SN'-') {MaTla + Ab76) 

= I ■ {{MaAab + kh)lh + MaBaP) . (112) 

As they are constraints, they have to vanish and consid- 
ering P and 7a as independent, it is possible to find Ba 
and Aab through 4 equations: 

EaBa = 0, MaBa = 0, (113) 

EaAab + Sb = 0, (114) 
MaAab+^b = (^- (115) 

With the expression of Q given in (jllOp . and using Eq. 
pi3p . one can see that, with our choice of generating 
function: 

= i?a{7a,H(2)(iveo)+p('Ha'^e)} 

= Ba{la,{N^o){EaTTa + ^blb)} 

+ Ba{7a,(a"C)(M,7r,+A,7,)}, 

= (TVCo) • BaEa + [d'^O ■ BaMa 

= 0. (116) 

This shows that Q2 is basically gauge-invariant because 
of the relations (|113p and not because the anomaly- 
freedom of the algebra. The gauge-invariance can also 
be seen by expressing Q2 in terms of {5E,Sip), and 
using (jlOOp . Finally, it is possible to define a set of 
3 new variables (Qb^Pb), and the function 5' by 5 = 
\Aablalb + BablaPb- Making this choice and applying 
the procedure described above leads to simple equations 
showing that Qb oc Q. This shows that there is a unique 
consistent choice for Q2 and the previous case is there- 
fore preferred. 



As far as the generating function is concerned, we could 
also have chosen S = fiija = q,Q) and found the con- 
jugate momentum P. This would however have led to 
Pi = —Q2 and the situation would have been equivalent. 
Moreover, it is also possible to consider a generating func- 
tion such that 

S = /3,4(P, {Q, P}) = ^AabTTaTTb + Ba^aP^^i. (117) 

Following the same procedure would lead to some new 
gauge-invariant functions. However they do not exhibit 
any interesting physical feature. In the following, we will 
therefore focus on the generating function given in (jlOOp , 
leading to Q = Q2 as a Mukhanov variable and P = P2 
as its momentum. 



C. Anomaly- freedom in the Hamilton- Jacobi 
approach 

To be consistent, that is to ensure that the evolution 
generated by the constraints remains compatible with the 
constraints themselves, the theory must be anomaly-free. 
In our case, using the same synthetic formulation as in 
piip and (|112l) . one can compute the Poisson brackets: 

{il(2)[^7Vi],ij(2)[^7V2]} = 0, 

{D'^^'^[5N^],d'-^^[SN^]} = 0. 

The total first order constraint M'^^^[SN, SN"] leads to: 

{A/(i)[l],Af(i)[2]} = [6N,6N-]{^aMa-EaAa). (118) 

To close algebra, that is to cancel anomalies, one has to 
require that 

J^aMa - EaAa = 0. (119) 

Using Eqs. PH)) and pT5)l : 

(HHl) X Mb = EaAabMb + ^aMa = 0, (120) 
dnSl) xEb = MaAabEb + KEa = 0. (121) 

Combining those equations with (jll9p . the condition for 
anomaly-freedom reads as: 

EaAabMb - MaAabEb, (122) 

which is fulfilled only if Aab is a symmetric matrix, with 
thus only 6 unknown parameters. This corresponds to 
a fully solvable problem. The condition of anomaly free- 
dom allows one to completely determine without ambigu- 
ity the equations of motion for the gauge-invariant vari- 
ables. 
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D. Mukhanov equation in General Relativity 



VI. SECOND CHANGE OF VARIABLES (Q,P) 



As perturbations can, a priori, be analyzed through 
different choices of gauges (for instance, the Newton 
Gauge where B = E = 0), it is useful to provide gauge- 
invariant quantities (related to 7a, 6N and 57V°) that 
are physically relevant to investigate observational con- 
sequences, the Bardeen Potentials p^ : 



In order to compute physical effects, one needs gauge- 
invariant variables. We have shown that the generat- 
ing function, defined in Eq. (|109|) . will give such gauge- 
invariant observables. In the following, we will precisely 
show that Q and P are the Mukhanov variables fulfilling 
the correct equations of motion. 



* = -r{B - E)+H{B - E), 



drj 

ip-H{B-E), 



(123) 
(124) 



where H is the conformal Hubble parameter. As we are 
dealing with a universe filled with a massive scalar field ip, 
it will also undergo gauge- invariant perturbations 6(p'~^^ . 
Gravity and matter perturbations are of course linked 
and wc shall focus on the linear order, as often in per- 
turbation theory. As derived in Q (this follows from the 
definition (|96)) ): 



5ip 



GI 



5ip + <^{B - E). 



(125) 



When dealing with all the scalar perturbations, there 
are 2 degrees of freedoms. This equation generates a 
constraint and only 1 degree of freedom remains. Usu- 
ally, the relevant variable used in cosmology is called the 
Mukhanov- Sasaki variable v, originally found in [ll| gov- 
erned by the associated Mukhanov-Sasaki equation: 



— + Aw - -w = 0, 
dri^ z 



where 



a{r]) 



a(ri) — . 



(126) 

(127) 
(128) 



When this variable is found, and after performing a 
Fourier transform, it is easy to compute the power spec- 
trum of the conserved curvature perturbation. As review, 
e.g., in [H: 



V = z ■ R, 



k^ 

Vnik) = - 



The spectral index is, for example, given by: 

d ln{PR) 



nc - 1 = 



d ln{k) 



(129) 



(130) 



(131) 



k=k^ 



This power spectrum PR{k, rf), typically representing the 
state of the universe at the end of infiation, is a manda- 
tory input to compute observables, in particular in the 
Cosmic Microwave Background (CMB). 



A. Expression of the gauge-invariant variables 



Using the requirement (jl07p with kc= 1 leads to con- 
ditions on Aab and Ba, through (|114p . They can be writ- 
ten as: 

H(') [SN] = = 1 • C + 70 • eo + 71 • 6 + 72 • 6, (132) 
where 





= 'J^B, 
a 






(133) 


eo 


= Aoo0 + 






(134) 


a 


= Aoi<^-|- 


a 




- (135) 
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= Ao2'f + 


2pk 

/I21 - 

a 


f - {2¥- - Kw) 


(136) 



Considering now Eq. (|108p . we are led to: 



where, 



^ + 70 • ^0 + 71 • ^1 4-72 ■ ^2, (137) 



(138) 
(139) 





= S2, 




So 


= b(p — 


A02 




ab k 


+ A 








K p 




S2 


= A22. 





(140) 
(141) 



This system is much simpler than in the ADM formal- 
ism and can be explicitly solved. Taking into account Eq. 
(jlOSp . one we can directly fix: 



Bo = A, 



0, 



An 



b(p, A12 



ab k 

K p 



(142) 



This choice for A02 and A12 leads to consider, in Eq. 
(|107p . only 3 equations for 5 unknown variables. One 
can check that the conditions (jl42p . implemented in Eq. 
P36p . make it vanish. This choice is therefore obviously 
correct. Moreover, as it will be made clear in the follow- 
ing, it is not necessary to determine all the coefficients 
of Aab- Let us now focus on Ba- Eq. (|133p gives one 
relation between the terms of Ba'. 



B-[ — — ^ - • Bp, . 
2pk " 



(143) 
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Bq will be kept as an irreducible degree of freedom and 
we will show that any gauge-invariant quantity will just 
be, at the end, proportional to the fundament Bq choice 



The new variable Q defined in Eq. (jllOp can be ex- 
pressed in terms of 7^ such that, finally, using 071 = 
— 2p'0, one obtains: 



= Bo ( (5</j - ^ • a7i 



= Bo[S^ 



(144) 



which is similar to v in (|127p if one chooses Bq = y/p = 
a{rf). It can also be noticed that it is independent of 
the choice of the base. Of course, one can also choose to 
invert Eq. (|TI51) . 



(145) 



2pk 

Bq — ^ • -Bi, 

and define the gauge-invariant variable as 

Q = Bi\ -70 + 71 = + i' 

\ a (fi J a \ip 

(146) 

which is a function of the perturbed curvature variable 
R: 



Q 



2pBi 



(147) 



which depends on a. In the following, we will focus on 
the case for which Q (x v but the other choice would also 
be possible and the same method would lead to similar 
results. Moreover, taking Bq = 1 instead of ^/p makes 
the calculation easier in order to find the function z, as 
in (|128p . In the next step, we will keep Bq free, until the 
last step, and derive the equations of motion for Q, and 
therefore for v, its Hamiltonian formulation and how to 
find the expression for z. 



Hamiltonian expression and Equations of 
Motions 



A general expression has been found for a gauge- 
invariant quantity Q which is related to the Mukhanov 
variables. As the generating functions and S are 
known, it is possible to find the Hamiltonian, and there- 
fore the Lagrangian, from which the equation of motion 
for Q can be derived. The classical results can then be 
obtained elegantly in the canonical formalism. Consider- 
ing 71 and 72 as pure gauge variables, as explained into 
the details in Q, we should avoid to use any function 
explicitly depending on them as they do not contribute 
to the dynamics. 

As we know the Hamiltonian as a function of {SK, 5E) 
and as we have derived the expression of the generat- 
ing functions, it is possible to express the second order 
constraints, that are governing the dynamics of perturba- 
tions, in terms of the new set of variables {Q,P). Using 



the notations of Q , the known variables can be inverted 
and one can easily find, with Bq = f{k,p): 



^'P = + [71,72], 

-DO 



A, 



<57r = i3oP+^Q + [71,72], 

TTl = —Q - —jBqP+ [71,72], 

Bq 2pk 



7r2 



'Bq 



[71,72] 



(148) 
(149) 
(150) 
(151) 



To go further in studying the dynamics, let us notice 
that, as Q is gauge-invariant. 



'JSN,SN' 



.Q^{Q,h''^\5N] + D^^^SN""]} ^Q. (152) 



The evolution of Q (and this is true for any gauge- 
invariant variable) is thus given by: 

Q = {g,iV(7j(") +-ff^^)) -H^"(£'(°^ +1?^^)) 

= + SsN,SN^Q + {Q,NH^^^ +N''D^^'>}. (153) 



So, 



Q = {Q,iVi/(2) +Ara£)(2)}. (154) 
In our case, these constraints are, with = 0: 

+ IvpS'^'daS^d.S^ 



(155) 



As SE is related to 7„i through Eq. (|41|) . we should not 
consider functions depending on it. Moreover, as 



— /(7m, 7m), 



(156) 



one can write, after taking the Fourier transformation 
and using Eq. ([7T|) : 



d^k^ 
2k 



3 ^1 



2 TTl TT2 



2 (be) 



+ +pV")-fl 



(157) 



The c parameter enters the equation only quadrati- 
cally, therefore the choice of the generating function S*^ 
(either 7 = P or 7 = Q) does not enter the final result. 
However, we do not have yet expressed the Hamiltonian 
for the gauge-invariant variables. It is necessary to use 
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Eqs. (|148p in (jl57p and to add the relevant terms as- 
sociated with the derivative with respect to time of the 
generating function, 

dS 1 • 



Bo 



To avoid inconsistencies, Q and P terms should not be 
taken into account. The gauge-invariant hamiltonian 
constraint for Q and P is thus: 



Hgi = 



2 



Bn 



p 




(159) 



where 



^00 

P 



An 



42 

^00 



3k 



-LP . 



(160) 



In order to recover the usual Hamiltonian formulation, 
the cross terms in P and Q in the previous equation 
should vanish. However, looking at (|159p . it is clear that 
this cross terms will give the expression of Aqq , as a func- 
tion of However, looking at F, we see that the square 
and the derivative with respect to time of Aqo are in- 
volved. For i?o 7^ 1, the expression is quite complicated 
(although in principle tractable) . As we are interested in 
the Mukhanov variable w, an easy choice is to first con- 
sider i?o = 1, to find the equation of motion for Q and 
its Hamiltonian, and then to set v ~ y/pQ and perform 
the related calculation. In the following, we will go on 
like this. Thus, to cancel the cross term, we set 



c^pk 



+ 2k- 



A 



00 



f> Aoo = - 



WithEq. (fmll : 



a 2k \ 2k 



(161) 



(162) 



Finally, with the expressions of Aqq, Aqi, and Bq, one 
can see that the Hamiltonian is independent of a and b, 
which means that the results is also independent of the 
normalization choice for the vectors A™. One can also 
solve (|135|) to find An but this is physically unuseful. 
One can also study the equations of evolution: once SE 
is known, it is naturally possible to obtain the exact 
value of SK by studying SE. This leads to valuable 
informations on tTq. 



Taking into account Eqs. (|16ip and (jl62l) . it is now 
possible to refine the expression of the Hamiltonian (jl59p 
as: 



Bgi 



d^k 



P 



(163) 



(158) 

where F is given by 



F = /s^ + py" + 3«;(^2 _^ intpV - - [ ) . (164) 

k 2 \ k ) 

To establish this expression, we have used the Raychaud- 
huri equation (written in conformal time iV = -^/p): 



n = k = k^^'-(p\ 

and the Klein-Gordon equation: 

(p^2k'^^-pV' = 0. 
The Hamilton equations thus lead to: 



Q 
p 



dHBj P 
dP p' 



dQ 



(165) 
(166) 

(167) 
(168) 



and we recover the usual Mukhanov-like equation: 

Q + 2kQ + rQ = 0, (169) 

remembering that k ~ Classically, the Mukhanov 
equation is derived from the action: 

1 



S= drj d-'k 



-k^ + - \ V 

z 



(170) 



The variable v was here defined with Bq = 1 such that 



leading, in the hamiltonian constraint, to: 



(171) 



_ ■ - P P 
V = kv + \ppQ = few H 7= — > — = = i) — few, (172) 

VP VP 



and so: 



B, 



GI 



d-^k 



liv-'kyf + ^Tv' 



(173) 



A standard Legendre transformation and an integration 
by part with respect to time leads to the Lagrangian: 



£^ = / d'k 



+ ^i-r + P + ky 



d-^k- 



• 2 I 2 2 I ^ 2 
V — k V H V 

z 



(174) 
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where terms leading to equivalent equations of motions 
have been deleted. With this equation, it is possible to 
recover z by solving; 



- = -r + 

z 



(175) 



which is satisfied classically for 



(176) 



This method allowed us to find gauge-invariant vari- 
ables and their equations of evolution, starting from an 
anomaly-free algebra in an easy way. In Loop Quan- 
tum Gravity, corrections to the classical theory are ex- 
pected and this will change the expression for the con- 
straints, leading to a modified algebra. The requirement 
of anomaly-freedom is not necessary to obtain the gauge- 
invariant Q and P, only a vanishing Poisson Bracket for 
the first order constraints is required. But, of course, full 
physical consistency can only be achieved if the algebra 
is closed. Thus, as in the previous sections where the 
case of General Relativity expressed with Ashtekar vari- 
ables was studied, we will consider in the following the 
effects of the two main corrections from LQG, that is the 
holonomy and the inverse- volume corrections. 



VII. APPLICATIONS 

We now consider constraints modified with respect 
to the classical case. In the following, we focus only on 
the steps useful to find the associated gauge-invariant 
Mukhanov variables, without going into the details of 
the calculations. The missing steps can easily be rebuilt 
using the techniques given above. 

We still consider a universe filled with a massive scalar 
field tp. The diffeomorphism constraint holds its classical 
form and, in all the following considerations, it will still 
be given by Eq. (|83p. Moreover, the expressions of 7^ 
and TTd do not rely on the shape of the constraints, but 
on the shape of the metric. In the following, expressions 
(|77ll80p will therefore be used. What will be modified 
are the Hamiltonian constraints where counter-terms 
have been added in order to cancel the anomalies so 
as to have a closed algebra. In the following, we will 
give the expressions of the first and second order for 
these constraints. The interested reader can go to 
the appendix where, starting from the zeroth order 
constraints, the equations of motion for the background 
variables are derived. 

To be as general as possible, we will keep the same 
notation where a and b are unknown. 



A. Inverse-volume case 

Following , where anomaly- freedom was found in the 
case of inverse- volume corrections, we consider hamilto- 
nian densities given by: 



1 



- — (l+g)kHEi- 



+ + 2^3/2 2p 



-3/2 



6^6 E'^ 

{l + f3)V^i^)S^ + V{^)^—^ 

2p 



In this case, the Friedmann equation is: 



— T9 i^y ly n _ 



and, by definition. 



3 V2p2 



(p = V—. 
P 



(177) 

(178) 
(179) 

(180) 
(181) 



Going ahead as in the classical case, Eq. (|108|) gives 
the relation between Bq and Bi such that: 



B jl + h) a ^ 
^ (1 + /) 2pak °' 



(182) 



Setting Bq = 1, the gauge- invariant Mukhanov-like vari- 
able is then: 



Q 



dS_ 
dP 



~ Sip - 



(1 + /i) 
(1 + /) ak 



(183) 



Proceeding as before to solve the Hamilton- Jacobi-like 
equations, Eq. (|136p has to be fulfilled, which can be 
expressed here as: 



6 



-aP(2/ + 5) + ^^(2/i 
2 V 



f2) 



In our approach, this condition is satisfied if 
9 = -2/, 

./2 = 2/1. 



(184) 



(185) 
(186) 



We have thus recovered exactly the relations given in 
so as to have an anomaly-free algebra. This is of 
course not surprising as (|184p is related to the condition 
of anomaly- freedom that was pointed out in (jll9p . 
In 0, the second order of the corrections a{p,SEf) and 
i^{p, SEf) also had to be taken into account, but as they 
are proportional to SE and thus to 71 and 72, we don't 
need to consider these terms. So, 



{a'^^\ 1^'^^^) = 7(71,72) ignored. 



(187) 
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The second order constraint density with N°- = can 
thus be written as: 

+ 2 + ^^^^ + 2 ^ ^ 55)v^fc2 V 

+ i(l +56)p'/V^^(^)'5^' + [<S£;,r]. (188) 

In this case, the cross-term PQ are vanishing if one im- 
poses the following condition on ^00 ^ 



An 



-K — 



y) 2k (1 + /) (l+5i)' 
The Hamiltonian can therefore be written as: 

2 



(189) 



p(i + gi)(^) +riviVpQ)^ 



(190) 



where 



IV ^ K 



-2a - + -a[-] +iy{l+gi)^ 



Let us define now 



D{l+gi) 



Q, 



vKi + ffi)- 



1 



The Lagrangian can be written as 

^GI 



— {v'+{-k'd{l+g^) + -]v^ 



(191) 

(192) 
(193) 

(194) 



where 



- = -Tiv + k^a{l + .95) + e + 



ak 



drj 



(195) 



Classically, the conserved curvature perturbation is given 
by: 



R^tp + -Sip = -. 

If z 



In the case of inverse- volume corrections, from the pre- 
vious equation and considering Eq. (jl92|) . one can natu- 
rally suggest in our approach: 



"(I + /1) 



P{l+gi) (1 + f) afc' 



(197) 



which is close but not exactly similar to the expression 
given in [l3l |. In fact, the propagation speed for the per- 
turbations derived in [l3| is given by: 



^ paper — (1 /s)) 

whereas, in our case, it is equal to: 

„2 -2 1- 



^'{1+95)- 



(198) 



(199) 



In this study, we have given some arguments to establish 
the function z. Although it would, in principle, be possi- 
ble to check its consistency, using Eq. (jl95p . this would 
lead to lengthy calculations that have not yet been car- 
ried out. It is however clear that our choice is associated 
with a correct Lagrangian. It might be that both solu- 
tions are physically equivalent. We let this question open 
for future studies. 



B. S7-LQC model — holonomy corrections 

We now focus on the case of holonomy corrections and 
we use the notations of [3] . 

The first order corrected constraints, with counter- 
terms ai not yet fixed but introduced to close the algebra, 
are given by: 



n 



(1) 



1 



__L(K [if + a2)5Ei + ^dcd'SE] ] , (200) 



p3/2 ~ 2p3/2 2p 



F53/2 



2p 



(201) 
(202) 



where we use the notation (n 7^ 0): 



K\ 



sin(n/i7fc) 



(203) 



One also has to deal with the Klein-Gordon equation. 



TT = pip, 

= -pd^V{^)-2K[2]0, 



(196) and the Raychaudhuri equation. 



fc = fcK [2] - f7 



0' 



with, 



ri= cos(2^7fc). 



(204) 
(205) 



(206) 



(207) 
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Moreover, holonomy corrections lead to a modified Fried- 
mann equation. The energy density p is basically defined 
through: 



and the Friedniann equation is given by: 



where 



H^ = IK[2]^ = -pp(^l-^), 



(208) 



(209) 



(210) 



is the critical energy density. Applying the same proce- 
dure as previously, one can derive the relation between 
Bq and Bi: 



Bi 



aip 



2p(K [si] + ai 



-Bo. 



(211) 



The related gauge-invariant variable is then (with the 
usual choice Bq = 1): 



(K [si]+ai) 
In this case, the condition p36p reads as 



(212) 



+2KnpV'- 



2K[2 

(p 1 / ^KCp 

C|2] ^ 2 VK[2] 



- = To - r! • fc' 

z 



\2X 



dK [2] 



which is fulfilled for 



"=K[2l^- 



(218) 



The Mukhanov equation in conformal time [N = y^), 
remembering that K [2] = is: 

Q + 2K [2] + TfiQ = 0. (219) 
As previously, it is possible to find z through 



(220) 



(221) 



This corresponds exactly to what was found, following 
another path based on the Bardeen potentials, in [l3 |. 
As in the previous case, we have recovered the Mukhanov 
Lagrangian. 



C. General case: inverse-volume and holonomy 
corrections 



= 2fc(IK[si] +ai) - 2IK[l]Va2, (213) 

which is again a necessary condition appearing when 
is considered and which has to be fulfilled to 
have an anomaly-free algebra. By this procedure, we 
have two unknown counter-terms a.\ and a2, and the 
previous equation gives a relation between them: 



a2 = 2K [1] - 2fc(IK [si] + a{). 



(214) 



In [iJli the anomaly was removed with a\ = K [2] — K [si] 
and thus = 2K [1] — 2fcK [2] which will be used in 
the following. The second order hamiltonian constraint 
density is thus modified such that: 



2k 



One obtains: 



P V 



2K[2y 



and: 



Am 
a 



2K [2] 



2K [2] J ' 



(215) 



(216) 



(217) 



In this section, we will not address the issue of 
anomaly-freedom for the case where both corrections 
are taken into account. We will just focus on defining 
the Mukhanov variable by the method previously devel- 
oped. Naturally, it then will be expressed as a function 
of counter-terms. We will see that this case can be solved 
as if corrections were mostly independent, as suggested 
in [3. 

In this case, the first order constraint densities can be 
defined, as in the previous case, such that: 



^G^ = ;^h4V^(K[,5i]+ai)5A-| 



VP 

for the gravitationnal sector and. 



, (222) 



SUE^ 



zp 



,(224) 



for the matter sector. 

The Friedmann-like equation is given by: 



(225) 
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with the definition: 



J. _7r 

'■P = • 
P 



(226) 



Going ahead as previously, one obtains again a relation 
between Bq and Bi. 



(1 + /i) a 



(1 + /) 2pa{K[si]+ai) 



Be, 



(227) 



As a fully closed algebra has not yet been found, we are 
compelled to keep ai as unknown. Using previous results 
from [l^ and jl4| . it is however probable that the solution 
will soon be derived. 

Setting i?o = 1, the gauge-invariant Mukhanov-like 
variable is then: 



Sip - 



(l + /i) 



(1 + /) a(K[si]+ai) 



ip. (228) 



Following the same procedure as in the previous sec- 
tions when solving the Hamilton- Jacobi-like equations, 
we have to fulfill, in particular, Eq. (jl36|) . which can be 
expressed here as: 



6 



6 



-a{2k{K [si] + ai) -t- - 2IK [1]^) (229) 



/2) 



and vanishes for 



2fc(K[si] +ai) + a2- 2IK[1]^ 
,/2 = 2/i. 



0, 
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(231) 
(232) 



This corresponds to conditions already found when 
holonomy or inverse- volume corrections were taken into 
account independently. 

At this stage, it is difficult to go much ahead. However, 
the previous results lead us to assume 



V 



i>{l+gi) (1 + /) a(K[si]+ai)" 
This remains to be fully demonstrated. 
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VIII. CONCLUSION 

This article builds on the innovative ideas given in 
and develop them so that they can be used in the 
framework of Loop Quantum Cosmology. Going through 
successive changes of variables by using appropriate 
generating functions, we have set a scheme useful for 
the definition of gauge-invariant variables appropriate to 
study cosmology. Both the cases of inverse- volume and 
of holonomy corrections were considered. In principle, it 
is also possible to study in such an easy way, any other 
kind of correction that can be set up by contraints. 



Moreover, the Hamilton- Jacobi method is very general 
and can be used in different frameworks, such as particle 
physics, where gauge-invariance is required. 

Several developments are expected. First, although 
reasonable, some guesses had to be made. They should 
be checked in details by going though the full exhaustive 
computation, in particular by studying the equations of 
motion for the Bardeen potentials. 

Second, the method should be applied again when 
counter-terms for the holonomy -f- inverse-volume case 
will have been found. This work is already on the way. 

Finally, those gauge-invariant variables should now be 
used to investigate cosmological consequences at the ob- 
servational level. 



IX. APPENDIX : EQUATIONS OF MOTION 

In this appendix, we give an example of how to de- 
rive the equations of motion for the background variables 
when holonomy and inverse- volume corrections are taken 
into account together. The classical limit corresponds to 
a — >■ or K [n] — )■ fc. At the lowest order, the diffeor- 
morphism constraints are null, we therefore consider in 
the following only the hamiltonian constraints which are 
given, for gravity and for matter, by: 



and 



— I d-'xN 
2k 



d-^xN V 



(-6aK [1]' Vl) , (234) 



2fV2 



The background variables are linked through: 



K 

3' 
1. 



(235) 

(236) 
(237) 



For instance, the conformal Hubble parameter found with 
p = a^irf) is given by: 



P_ 
2p 




dk 

(-12V?aIK [2]) 



aK [2] 



Moreover, the energy density p is defined by 



1 5H„ 



p-2 



6N 



2 p3 



V, 



(238) 

(239) 

(240) 
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(242) 
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and is linked to the gravity through an equation of motion 

(243) 

which gives the Friedniann-hke equation 



aK [ly 



2 p2 



pV 



(244) 



Using the previous relations, the Fricdmann equation is 
thus given by: 



{aK [2])2 = a{aK [if) ■ cos^{fijk) (245) 



«(aK[lf)(l-(/27)']K[lf) 
V 3 / \ a pc 



K- / - P 

-PP " 

P 



(246) 
(247) 

(248) 



where, as usual when using the holonomies within the 
/i-scheme, the critical energy density is defined by Eq. 
(|210p . Moreover, similarly, the Raychaudhuri equation is 
given by 



akK [2] 



^SaK[l]^ 
a op / 



K If 

2 V 



(249) 
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